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In this paper we examine the effects of electron-hole asymmetry as a consequence of strong cor- 
relations on the electronic Raman scattering in the normal state of copper oxide high temperature 
superconductors. Using determinant quantum Monte Carlo simulations of the single-band Hub- 
bard model, we construct the electronic Raman response from single particle Green's functions 
and explore the differences in the spectra for electron and hole doping away from half filling. The 
theoretical results are compared to new and existing Raman scattering experiments on hole-doped 
La2- x Sr x Cu04 and electron-doped Nd2- x Ce x CuC>4. These findings suggest that the Hubbard model 
with fixed interaction strength qualitatively captures the doping and temperature dependence of the 
Raman spectra for both electron and hole doped systems, indicating that the Hubbard parameter 
U does not need to be doping dependent to capture the essence of this asymmetry. 

PACS numbers: 78.30.-j, 74.72.-h, 71.10.Fd, 74.25. nd 



I. INTRODUCTION 

The parent compounds of cuprate high temperature 
superconductors are antiferromagnetic Mott insulators at 
half filling. 1 By removing electrons from or adding elec- 
trons to the CuC>2 planes by chemical substitution the 
antifcrromagnctism is suppressed and superconductivity 
appears over a limited doping range. 2 At first glance 
one would expect that the doping leads to effects more 
or less symmetric around half filling similar to the re- 
cently discovered FeAs superconductors. 3 However, in 
the cuprates the differences originating from cither elec- 
tron or hole doping can be quite significant: the max- 
imal superconducting transition temperature T c hardly 
exceeds 30 K for electron-doped cuprates while reach- 
ing 150 K for hole-doped materials; in the normal state, 
while the approximately linear variation with T of the 
resistivity over wide temperature ranges is a hallmark of 
the hole-doped systems there is much more doping de- 
pendence on the electron-doped side where the resistiv- 
ity crosses over to a nearly T 2 behavior already slightly 
above optimal doping close to x = 0.15; in the Raman 
spectra of hole-doped systems, the B^g response is essen- 
tially universal over the entire doping range of the su- 
perconducting dome whereas, concomitant with the re- 
sistivity, the Raman spectra of electron-doped systems 
changes rapidly and, at low temperatures, exhibits Fermi 
liquid-like shapes for x > 0.16. 

Nevertheless, this asymmetry is not entirely unex- 
pected for strongly correlated copper oxides. 1,11,12 In 
the insulator at half filling, the wavefunction is com- 
posed largely of a superposition of a Cu d 9 hole on 



each 3d x 2_ y 2 orbital, with a minority of d 10 L charac- 
ter, wherein the Cu orbitals are filled and a hole occu- 
pies the oxygen 2p ligand L. When hole-doped away 
from half filling, a d 9 L state forms wherein the addi- 
tional hole gains derealization energy as well as mag- 
netic exchange energy by occupying the oxygen ligand 
to form a so-called Zhang-Rice singlet. 13,14 In contrast 
doped electrons tend to reside solely on Cu d 10 sites. This 
asymmetry can be revealed by comparing photoemission 
with inverse photoemission, or via angle-resolved pho- 
toemission spectroscopy (ARPES) in hole- and electron- 
doped cuprates. 10,15 More recently the issue of particle- 
hole asymmetry has been well explored in scanning tun- 
neling microscopy (STM) studies. 16-19 

In this paper we explore how particle-hole asymme- 
try can be viewed from Raman scattering measurements. 
In particular, we are motivated to explore the ques- 
tion of whether the low energy particle-hole excitations 
emerge from doping a Mott insulator while preserving 
the strength of correlations in the undoped parent in the 
form of Hubbard U, or if these excitations are better 
described in terms of a strongly doping-dependent Hub- 
bard U leading naturally to a collapse of the Mott gap 
not driven by simply adding particles, but by a strong 
decrease in U with doping. Recent comparisons of spec- 
tral weight transfer observed with ARPES, x-ray absorp- 
tion, and optical spectra have been interpreted in terms 
of a doping dependent U. 20 This differs from the con- 
clusion reached from many other studies of the Hubbard 
model. 21 - 24 

We construct the electronic Raman response using 
single-particle propagators determined from determinant 
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quantum Monte Carlo simulations of the single-band 
Hubbard model. The theoretical results that highlight 
differences in the spectra for hole and electron doping 
are compared to results of Raman scattering experiments 
on La2- x Sr x Cu04 and Nd2- x Ce x Cu04. This qualitative 
comparison suggests that a constant Hubbard interaction 
U captures the essence of the doping and temperature 
dependence of the particle-hole asymmetry. 

In Section II we present a brief description of the theo- 
retical calculation and the main results. Section III pro- 
vides details on the sample preparation, the experimental 
methods and results for comparison to theoretical cal- 
culations. A discussion of qualitative similarities and 
differences between experiment and theory appears in 
Section IV including comparisons between the extracted 
scattering rates (Raman resistivity) and the evolution of 
the Raman spectral weight with doping and temperature. 
Finally, we present conclusions in Section V. 



II. THEORY 
A. Model 

The single-band Hubbard Hamiltonian represents an 
effective low-energy model for the cuprates. 13 ' 25 Its ap- 
plicability derives from down-folding models explicitly 
incorporating planar copper and oxygen degrees of free- 
dom. Written in a second-quantized real-space represen- 
tation, the Hamiltonian takes the form 

lj,<7 i,(7 

+ ^EKt-^)(«a^)- (i) 

i 

The operators c] a and Ci ](T create or annihilate an elec- 
tron with spin a at site i, respectively, and rii tC , = c\ a Ci t(T 
in each spin channel. The non-zero tight-binding coeffi- 
cients {Uj}, restricted to nearest-neighbor t and next- 
nearest-neighbor t' hopping, together with the chemical 
potential fi, that controls the electron filling, define the 
nonintcracting bandstructure and the Hubbard repulsion 
U controls the strength of electron-electron correlations. 

While this Hamiltonian appears rather simple, it re- 
sists an analytical solution in two-dimensions, applica- 
ble to the cuprates, and is challenging to solve numer- 
ically, especially for the intermediate range of interac- 
tion strengths U believed to be most appropriate to the 
cupratc problem. We choose to work with U = 8t, 
equal to the noninteracting bandwidth W, that repre- 
sents a canonical value for the interaction strength in 
the cuprates related to the charge-transfer energy be- 
tween copper and oxygen orbitals in these systems; it 
also sets the largest energy scale for the problem that 
can be observed directly in the high energy Raman re- 
sponse. Throughout the theoretical analysis, the nearest- 



neighbor hopping t serves as the energy unit of the prob- 
lem and we substitute a reasonable estimate for this 
down-folded hopping integral only for the purpose of 
comparison to experimental results. 

We numerically investigate the Hamiltonian of Eq. (I) 
using determinant quantum Monte Carlo (DQMC), 26 ' 27 
an auxiliary-field technique. In principle, this method 
is numerically exact and allows one to determine both 
single- and multi-particle response functions at finite 
temperature. However, computational costs limit inves- 
tigations to finite-size, small clusters (in either real- or 
momentum-space) and the fermion sign problem 28,29 lim- 
its the lowest accessible temperatures where one can still 
obtain reasonably accurate results. The finite-size clus- 
ters used in this study are larger than those that can be 
accessed using exact diagonalization and provide a suf- 
ficient sampling for reconstructing details of the single- 
particle self-energy assumed to vary slowly as a function 
of momentum for the chosen parameters. 

The DQMC method supplies the finite temperature, 
imaginary time propagator Gij (r) on a finite-size cluster 
with periodic boundary conditions. Individual Markov 
chains of the Monte Carlo process provide input for the 
maximum entropy method (MEM) 30,31 used to Wick ro- 
tate the imaginary time data to real frequencies using 
Bayesian inference from separate estimates of the prop- 
agator assumed to have a Gaussian statistical distribu- 
tion over different chains. The data are characterized 
to ensure that they reasonably satisfy this assumption; 
and they are preprocessed and/or more data are gath- 
ered to satisfy these conditions. From the real-space 
statistical ensemble {Gjj(r)}, a discrete Fourier trans- 
form yields {Gk(t)} from which MEM returns the single- 
particle spectral function A(JZ, uj) on the corresponding 
discrete momentum grid. Once obtained in this fashion, 
the single-particle self-energy E(K, uj) can be extracted 
using Dyson's equation and the bare bandstructure corre- 
sponding to the tight-binding model parameters. Assum- 
ing a weak momentum dependence to the self-energy, an 
interpolation routine provides the value of E(k, w) at an 
arbitrary point k in the Brillouin zone (BZ) and Dyson's 
equation can be employed to compute A(k 7 uj) at that 
point. 

Two-particle response functions such as the charge 
and spin susceptibility (the dynamic structure factors 
S(q, Q) or optical conductivity er(Q)) can be evaluated 
in imaginary time and analytically continued to real fre- 
quency using a similar prescription to that followed for 
the single-particle Green's function. 30 These quantities 
satisfy well-defined sum rules that make redefining the 
spectral functions in terms of probability distributions 
and subsequently normalizing the imaginary time data 
relatively straightforward. While in principle bounded, 
the Raman response does not satisfy any similar sum 
rule; and considering the significant fermion sign prob- 
lem that already complicates the analytic continuation 
by adding an additional source of noise and covariance in 
the data, rather than evaluate the Raman response for 
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different symmetries directly in imaginary time (or Mat- 
subara frequency) , we evaluate the single-particle Green's 
function in imaginary time, Wick rotate using MEM, and 
then estimate the Raman response as 8 

oo 

= T^E^W / G"(k, w )G"(k,c + 0) 

K -oo 

x [/(w) - f(u> + fi)] duo. (2) 

In practice the integral over real frequencies is evaluated 
numerically using Ricmannian integration and the up- 
per and lower limits of the integral are cut-off at fre- 
quencies > 5t beyond the "step-edges" set by the dif- 
ference in Fermi functions appearing in the integrand, 
providing sufficient accuracy over the studied tempera- 
ture interval. The vertices 7/j(k) are chosen to corre- 
spond to Big [7s lg (k) = \(cos(k x ) - cos(fcy))] and B 2g 
[7 B2s (k) = sin(A;a;)siri(fcj / )] symmetries that highlight the 
anti-nodal and nodal portions of the Fermi surface, re- 
spectively. While this method neglects vertex correc- 
tions in the Raman response, it captures features that 
correlate with different intra- and interband charge ex- 
citations within the model that qualitatively compare 
to results from experiments on the cuprates. The rela- 
tively high temperatures, lack of vertex corrections, and 
the simplified Hamiltonian mean that other low energy 
features that can be seen in the experiment like multi- 
magnon excitations and phonon degrees of freedom are 
missing from this analysis. While the appearance of 
multi-magnon excitations in the Raman response is usu- 
ally attributed to the effects of higher-order resonant dia- 
grams off-resonance, 32,33 even two-particle vertex renor- 
malization should contribute to the appearance of these 
features in the response. 34 

We use 64-site square clusters with periodic boundary 
conditions corresponding to a momentum space grid {A"} 
with spacing 7r/4 in each direction. The imaginary time 
interval has been partitioned into L = 48 "slices" of size 
At = /3/L running from to /3. As noted t serves as 
the energy unit of the problem. For this study, j3 varies 
between 1/t and 3/t giving a value of Ar that varies 
between 1/48/; and 1/164, controlling the Trotter error 
while maintaining a reasonable computational time to 
completion, with the majority of results shown for /3 = 
3/t. 

B. Results 

Before exploring the Raman response for different re- 
gions of parameter space for the Hamiltonian of Eq. (1), 
let us first look at the single-particle spectral function for 
the half-filled model to understand the nature of inter- 
band charge excitations that can appear in the response 
function. Figure 1 displays the calculated bandstructure 
of the half-filled single-band Hubbard model with param- 
eters t' = -0.30<, [i = O.OOi, fj = 3.0/t along high sym- 
metry directions in the BZ. Immediately noticeable are 
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FIG. 1: (Color online.) Theoretical band dispersion along 
high symmetry directions for the single-band Hubbard model 
at half-filling with parameters t' = — 0.3/;, fi = 0.0/;, U = 8.0t, 
and P = 3.0/t obtained using determinant quantum Monte 
Carlo as described in the main text. The red (solid) arrows 
highlight possible q = (0, 0) transitions in the Bi g Raman 
scattering channel while the purple (dashed) arrows highlight 
possible transitions in the i?2g channel along these high sym- 
metry cuts. Adapted from Ref. 35. 



the incoherent lower and upper Hubbard bands (LHB 
and UHB) centered near the T-point and (7r,7r), respec- 
tively. Above the Fermi level, the UHB has a dispersing 
branch along the (tt, 0) — (n, tt) and (0, 0) — (tt, tt) direc- 
tions. Along the (0, 0) — (tt, 0) direction, this feature is 
nearly dispersionless. In particular, there is significant 
spectral weight in the region near (tt, 0) at binding en- 
ergies near 2 t. Below the Fermi level, while the bulk 
of the LHB spectral weight is concentrated near the T- 
point, there is a dispersing branch along the (0, 0) — (7r, 0) 
and (0, 0) — (7r, 7t) directions that appears to be most pro- 
nounced at binding energies near —It. These features are 
qualitatively similar to those observed in experiment 36 
where a dispersive feature near (tt/2,tt/2) crosses-over 
to the higher energy valence band, assumed to have sig- 
nificant oxygen character. These dispersing features in 
the LHB and UHB are precursors to a quasiparticle-like 
band crossing the Fermi level that appears upon either 
hole or electron doping. 

The vertical double-headed arrows (red (solid) in the 
Big channel and purple (dashed) in the B 2g channel) that 
appear in Fig. 1 mark the energy scale of possible inter- 
band charge excitations that can be observed in the var- 
ious Raman scattering channels. Note that the Raman 
Big and B 2g vertices highlight the anti-nodal and nodal 
regions of the BZ, respectively, and, by symmetry, are 
identically zero along certain high symmetry directions 
as revealed in the form of each vertex entering Eq. (2). 
The lowest energy scale for each symmetry is associated 
with the insulating Mott gap with an onset energy ~ 2 1 
and a weak tail at lower energies due to the relatively 
high temperature. The main peak associated with the 
interband transition across the Mott gap should occur in 
both channels at an energy ~ 4 1 with a further promi- 
nent transition between the LHB and dispersing tail of 
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FIG. 2: (Color online). Theoretical band dispersion along 
high symmetry directions for the single-band Hubbard model 
near optimal (~ 15%) hole-doping (t' = — 0.3t, /j, = — 2.5t, 
U = 8.0t, and j3 = 3.0/t) with red (solid) and purple (dashed) 
arrows highlighting the prominent transitions in the Raman 
response in Bi g and B29 channels, respectively. 
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FIG. 3: (Color online). Theoretical band dispersion along 
high symmetry directions for the single-band Hubbard model 
near optimal (~ 15%) electron doping (t 1 = — 0.3f, fi = 2.0t, 
U — 8.0t, and /3 = 3.0/t) with red (solid) and purple (dashed) 
arrows highlighting the prominent transitions in the Bi g and 
B2 g Raman response, respectively. 

the UHB (or precursor to the quasiparticle-like band) at 
energies between ~ 6t and ~ 10 t, although this could be 
fairly broad. 

Upon either hole or electron doping, the gap at the 
Fermi level closes and the chemical potential moves into 
the dispersive portions of either the LHB or UHB forming 
a quasiparticle-like band (QPB) at low binding energies 
near the Fermi level. This would be reflected in the Ra- 
man response by an onset directly at zero energy and 
significant low energy (quasiparticle) spectral weight. 

Figure 2 shows the calculated bandstructure for the 
hole-doped single-band Hubbard model near optimal 
doping (~ 15%) with parameters t' = — 0.3t, /z = — 2.5t, 
and /3 = 3.0/t, again along high symmetry directions 
in the BZ. The QPB here is obviously well separated 
from the LHB given the significant coexistence of both 
features along the (0,0) — (tt, 7t) and, even more so, the 
(0, 0) — (tt, 0) directions. The QPB crosses the Fermi level 



near (7r/2, 7t/2) and at ~ (vr/4, 7r/4) the spectral intensity 
drops, demarcating a cross-over between the QPB and 
the LHB. On the whole, the evolution of the QPB qual- 
itatively agrees with the results of ARPES experiments 
on hole-doped compounds, 37-41 including the evolution 
of spectral intensity and changes in momentum space po- 
sition and robustness of this "waterfall" -like appearance 
as a function of momentum as found in previous work on 
the single-particle bandstructure. 22,35,42-46 

Figure 3 shows the spectral function for an electron- 
doped system with model parameters t' = — 0.3i, fi = 
2.0t, and fi = 3.0/t near optimal electron-doping (~ 
15%). The LHB, centered at ~ — 6i, has been reduced 
in intensity from spectral weight transfer into the QPB 
which now disperses down across the Fermi level from the 
precursor in the UHB. The QPB reaches approximately 
twice as far below the Fermi level than the QPB under 
hole-doping. This dichotomy or asymmetry in the band 
dispersion between hole and electron-doped systems can 
be traced back to differences in the shift of the chemical 
potential either into the LHB or the UHB with doping 
and the character of the state that then disperses across 
the Fermi level. 35,46 

In the hole-doped system the lowest energy scale de- 
tected in the Raman response should reflect intraband 
transitions within the QPB close to the Fermi level at an 
energy near 1 1. This energy scale also partially reflects 
the fairly high temperature of this study and would likely 
decrease with reduction in the model temperature. An 
interband transition between the LHB and QPB occurs 
between ~ 2 1 and 3 1 with yet a higher energy transition 
possible between the LHB and UHB at energies between 
~ 6 1 and 9 1. The features in the Raman response should 
all be fairly broad, not only because of the high simula- 
tion temperature, but also because of the intrinsically 
broad incoherent LHB and UHB that has a diminished 
spectral weight due to transfers into the QPB. The tran- 
sitions that should be prominent in the response are indi- 
cated by red (solid) and purple (dashed) arrows in Fig. 2 
for the Bi g and B-i g scattering channels, respectively. 

In the electron-doped system the general reduction in 
the LHB spectral weight and the lack of significant dis- 
persive portion of the UHB at higher energies above the 
Fermi level reduce the number of prominent transitions. 
The double-headed arrows in Fig. 3 indicate two promi- 
nent transitions: one from intraband transitions at low 
energies under 1 1 and an additional one from interband 
transitions between the weak incoherent LHB and dis- 
persive QPB near It. A weak tail exists in the UHB at 
higher energies; therefore, interband transitions from the 
QPB or LHB to the UHB also would be weak. 

Figure 4 shows the response calculated using Eq. (2) 
in both the B\ g and B^ g scattering channels for the half- 
filled single-band Hubbard model. Both channels show 
an onset at ~ 2 1 corresponding to the Mott gap present 
in the half-filled model with a weak tail at lower energy, 
due to finite temperature effects, that should shrink at 
lower temperatures. Between 4 1 and 5 1 there is a promi- 
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FIG. 4: (Color online). Theoretical Raman response in B\ g 
(red, solid curve) and £?2 9 (purple, dashed curve) symmetries 
for the half-filled single-band Hubbard model. (Parameters 
as in Fig. 1). 
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FIG. 5: (Color online). Theoretical Raman response in the 
[(a) and (b)] B\ g and [(c) and (d)] B2 g channels for various 
doping levels on the [(a) and (c)] electron and [(b) and (d)] 
hole sides of the phase diagram for the single-band Hubbard 
model. 



nent peak in the Bi g response and a shoulder in the B\ g 
response corresponding to transitions between the dis- 
persing portions of the LHB and UHB. This gives way to 
the strong peak in B\ g symmetry and broad shoulder in 
Big symmetry at approximately 8 1, associated with ad- 
ditional intcrband scattering pathways as indicated pre- 
viously. 

Figure 5 shows the B\ g and Bi g Raman response for 
various values of electron [panels (a) and (c)] and hole 
[panels (b) and (d)] doping. Both B\ g and B^ g symme- 
tries show progressive transfer of spectral weight to lower 
energies with doping away from half-filling as the intra- 
band transition dominates and the single-particle spec- 
tral weight is transfcred to the QPB. At low hole-doping 
in = 0.95) the intraband transition appears as a small 
shoulder or knee at low energy and the B\ g and Bi g re- 
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FIG. 6: (Color online). Theoretical low-energy Raman re- 
sponse in [(a) and (b)] Bi g and [(c) and (d)] B^ g symmetry 
from Fig. 5. These results have been converted to cm -1 as- 
suming a reasonable value of t = 400 meV that provides a 
qualitatively good description of the single-particle spectral 
function in these systems. 



sponse arc dominated by the high energy transitions from 
the LHB to the UHB. Near optimal doping (n = 0.85) the 
Raman response behaves as indicated previously with the 
low energy response in the Bi g channel dominated by the 
LHB to QPB transition that also gives a very broad peak 
in Bi g symmetry. With additional overdoping (n = 0.75) 
the response is dominated by the intraband transition 
with the LHB to QPB transition giving an asymmetric 
shoulder in B± g and large broad peak in Big symmetry. 
Upon electron doping, the Raman response in each chan- 
nel is dominated by the intraband transition at low en- 
ergy and the LHB to QPB transition at higher energies. 
While the strength of each feature varies with doping, the 
energy position remains relatively fixed with increasing 
electron count. 

Figure 6 shows the Raman response in both channels at 
low energy associated primarily with the intraband tran- 
sition. The energy scale has been expressed in cm -1 as- 
suming a value of t = 400 meV. For hole doping, the low 
energy peak grows and appears to shift in energy from 
- 2000 cm" 1 at low doping (n = 0.95) to - 4000 - 5000 
cm -1 on the overdoped side of the hole-doping phase di- 
agram (n = 0.75). However, this apparent shift of energy 
scale is presumably due to the overlap of several peaks 
that possibly could be distinguished upon lowering the 
temperature. With electron doping, the peak intensity 
grows and the energy shifts to slightly smaller values at 
the largest electron count (n = 1.20) near ~ 2000 cm -1 . 
In both cases, upon reducing the temperature from that 
used in the simulation one expects that the peak will 
narrow and shift to even lower energies. 

The low energy Raman response, particularly as Q — > 
0, can be used to obtain an estimate for the effective 
quasiparticle scattering rate in each channel and high- 
light important differences between hole- and electron- 
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doped systems. This Raman resistivity can be deter- 
mined from 8 



r M (T) 



dVL 



-i -l 



lim 



n 



(3) 



Figures 9(c) and (d) show the Raman resistivity in the 
Bi g and B 2g channels for both hole- and electron-doped 
systems near optimal doping extracted from the data pre- 
sented in Fig. 6. As one may expect, the effective scat- 
tering rate decreases with decreasing temperature as well 
as increasing doping away from half- filling (not shown) 
on either the hole- or electron-doped sides of the phase 
diagram indicating that the doped system becomes pro- 
gressively more metallic. From Fig. 6 one also can see 
that the slope of the Raman response is always larger for 
the electron-doped models compared to their counter- 
parts with similar hole doping. This behavior is further 
reflected in the Raman resistivities shown in Figs. 9(c) 
and (d). 



III. EXPERIMENT 

A. Samples and experimental details 

Experiments were performed on single crystals of 
La 2 - x Sr x Cu0 4 (LSCO) and Nd 2 - x Ce x Cu0 4 (NCCO) on 
the hole- and electron-doped sides of the phase diagram, 
respectively. The doping level x is indicated in each fig- 
ure. Both compounds belong to the 214-family and crys- 
tallize in the T'- and T-structure with and without oxy- 
gen atoms in the apex position, respectively. The crys- 
tals were prepared via the traveling solvent floating zone 
(TSFZ) technique. The Nd 2 - x Ce x Cu0 4 samples had to 
be post-annealed in pure Argon to remove the excess oxy- 
gen at the apex position and make the samples supercon- 
ducting. 47 ' 48 The annealing protocols for La2- x Sr x Cu04 
are described elsewhere. 48 



The Raman experiments were performed with a stan- 
dard light scattering setup. For excitation an Ar + laser 
was used and operated at 458 nm for La2- x Sr x CuC>4 and 
514 nm for Nd2- x Ce x Cu04. The angle of incidence of 
the exciting light was 66° . To achieve proper polarization 
states inside the sample the polarization of the light out- 
side was controlled with a Solcil-Babinct compensator. 
The samples were mounted on the cold finger of a He- 
flow cryostat with temperatures in the range from 4 to 
330 K and a vacuum of better than 10 -6 mbar. The scat- 
tered light was collected with an objective lens. Photons 
with selected polarization states were analyzed using a 
Jarrell-Ash 25-100 scanning spectrometer equipped with 
a CCD camera. The resolution at 458 nm was 9.5 cm -1 
unless otherwise stated. All spectra are divided by the 



thermal Bosc factor {1 + n(0,T)} = (1 



-n/T\ 



and 



corrected for the sensitivity of the entire setup including 
the energy dependence of the spectral resolution. 

By properly selecting the polarizations of the incident 
and scattered photons excitations of specific symmetries 
can be projected out. For particle-hole excitations in the 
cuprates, the B\ g spectra project mainly the principal 
axes while the B 2g spectra contain information about the 
diagonals of the tetragonal Brillouin zone. 50 The respec- 
tive light polarizations and Raman vertices are indicated 
as insets in Fig. 7. 



B. Results 

Fig. 7 shows the experimental B\ g and B 2g Raman 
spectra of Nd2- x Ce x Cu04 and La2- x Sr x Cu04 for ener- 
gies up to 5600 cm" 1 at 200 K. For undoped NCCO 
(x = 0.00), the most prominent peak in B\ g symmetry 
[Fig. 7(a)] is observed at 2900 cm -1 and originates from 
nearest-neighbor spin flip excitations 32, 49 > 51 ~ 53 which are 
not part of the theoretical description in Section II. At 
1200 cm" 1 there is a weak band that originates from two- 
phonon scattering. 54 Since the resolution is 28 cm -1 close 
to the laserlinc at 514 nm most of the phonon excita- 
tions are hardly visible. With doping the two-magnon 
excitation is suppressed rapidly but traces thereof may 
still be present at x = 0.12 (see also Ref. 55). Below 
2000 cm -1 there is little change of the continuum with 
doping. At high energies the intensities of the spectra do 
not depend in a systematic way on doping in contrast to 
what is observed for YBa2Cu306+a; (Y-123) 56 and also 
in LSCO as shown below. We assume a strong contribu- 
tion from luminescence, which can be seen by comparing 
annealed and as grown samples (for a discussion see e.g. 
Ref. 49), due to charge traps in the rather imperfect Nd- 
Ce-0 layers which mask the intrinsic effects of the carrier 
dynamics. 

In B 2g symmetry [Fig. 7(b)] the response below ap- 
proximately 2000 cm" 1 is weak for the undoped com- 
pound presumably due to the small carrier concentra- 
tion. The peak at 2900 cm -1 does not originate from 
polarization leakage otherwise the two-phonon excitation 
at 1200 cm -1 present in B\ g symmetry would be strong 
enough to be visible as well. The energy of the B 2g fea- 
ture with respect to that in B\ g symmetry is smaller than 
in LSCO. Since the next- nearest neighbor coupling J' de- 
termines by and large the peak energy in B 2g symmetry 52 
we conclude it is weaker in NCCO than in hole-doped sys- 
tems. The two-magnon scattering in B 2g symmetry dis- 
appears faster with doping than in B\ g symmetry. With 
increasing doping level the spectra gain intensity in the 
low as well as in the high energy range without signifi- 
cantly changing the spectral shape. 

In contrast, a strong doping dependence of the spec- 
tral shape is found for LSCO as shown in Fig. 7(c) and 
(d). The pronounced B\ g peak at 3300 cm" 1 observed 
at x — corresponds to two-magnon scattering and 
is progressively suppressed upon doping. The peak at 
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TABLE I: Complete list of samples studied partially adapted from Ref. 49. The results on La2- : rSr a ;Cu04 (LSCO) have been 
published in Ref. 49. Those on Nd2- :E Ce :!; Cu04 (NCCO) were taken on freshly prepared single crystals. 48 In the case results 
similar to ours were published before we give the references in the text and in the figure captions. Samples labeled with a have 
been prepared by M. Lambacher and A. Erb (WMI Garching), 48 b by Seiki Komiya and Yoichi Ando (CRIEPI, Tokyo and 
Osaka University), and c by N. Kikugawa and T. Fujita (Hiroshima and Tokyo). The transition temperatures were measured 
either resistively or via magnetometry or via the non-linear ac response. The T c of the 5% sample is the onset point of the 
transition. Tm was not measured for La2- x Sr x CuC>4 at x — 0.02 and 0.05. In the latter case Tjv = 0. 



c a rri nip 


sample ID doping 


T /TVr 
± c 1 ± N 


(K) AT (K) 


rrtTTiTTiPTlf 




La 2 Cu0 4 


LCO-00 


0.00 


0/325 


- 


Ar annealed 


a 


Lai.98Sro.o2CuC>4 


LSCO-02 


0.02 


0/- 


- 


as-grown 


<: 


Lai.gsSro.osCuCU 


LSCO-05 


0.05 


5/0 


3 


O2 annealed 


a 


Lai.92Sro.o8CuC>4 


LSCO-08 


0.08 


18 


4 


O2 annealed 


c 


Lai.ssSro.isCuO^ 


LSCO-15 


0.15 


38 


3 


O2 annealed 


a 


Lai.83Sro.i7CuC>4 


LSCO-17 


0.17 


39 


1 


O2 annealed 


b 


Lai.8oSro.2oCuC>4 


LSCO-20 


0.20 


24 


3 


as-grown 


a 


Lai.75Sro.26CuC>4 


LSCO-26 


0.26 


12 


3 


O2 annealed 


c 


Nd 2 Cu0 4 


NCCO-00 


0.00 







Ar annealed 


a 


Ndi.88Ceo.i2Cu04 


NCCO-12 


0.12 







Ar annealed 


a 


Ndi.s7Ceo.i3Cu04 


NCCO-13 


0.13 


9.9 


7.5 


Ar annealed 


a 


Ndi.8sCe .i5CuO 4 


NCCO-15 


0.15 


23.6 


1.3 


Ar annealed 


a 


Ndi.84Ceo.ieCu04 


NCCO-16 


0.16 


16.3 


2.5 


Ar annealed 


a 


Ndi.83Ceo.i7Cu04 


NCCO-17 


0.17 


5.0 


3.5 


Ar annealed 


a 



4200 cm" 1 which appears in B\ g and Bi g symmetry with 
the same intensity indicating its A^ g nature comes from 
higher order spin excitations including cyclic exchange of 
spins. 57 The low energy response is weak but picks-up in- 
tensity with doping. Here, it appears as if spectral weight 
would be transferred from high to low energies. In the 
high energy part the intensity increases between x = 
and 0.02 then decreases monotonically with doping. This 
results in a fairly fiat spectrum at high energy for the 
highest doping levels. Above x = 0.2 there appears a 
peak at low energy which is related to long-lived particle- 
hole excitations. This transfer of intensity is qualitatively 
predicted already in the Falicov-Kimball and Hubbard 
models. 58-61 A discussion in terms of a Fermi-liquid ap- 
proach is presented elsewhere. 62 

In Big symmetry, there are pronounced phonon bands 
below 1000 cm -1 in the undoped compound [Fig. 7(d)] 
which disappear quickly with doping. In contrast to the 
Big channel there is no redistribution of spectral weight 
from high to low energies upon doping. The overall vari- 
ation of the intensity is also non-monotonic. For the 
highest doping levels the continuum hardly depends on 
energy. 

It is particularly instructive to compare the temper- 
ature dependence at low energies of electron- and hole- 
doped materials. The Bi g Raman response of overdoped 
NCCO and LSCO are plotted in Fig. 8(a) and (b). In 
either case the initial slope increases upon cooling cor- 
responding to metallic behavior, however, with a signif- 
icant difference. In NCCO (x = 0.17) an isolated peak 



appears at low energy and temperature accompanied by 
a suppression of spectral weight in the range between 200 
and 800 cm -1 as already observed earlier by Koitzsch et 
al. 63 In overdoped LSCO (x = 0.26) there is only an over- 
all reduction of intensity in the entire range without any 
pile- up at low energies (see also Ref. 49). At tempera- 
tures above 200 K the spectral shapes become similar on 
both sides of half filling. 



IV. DISCUSSION 
A. Overall Features 

It is clear from the theoretical results presented in Sec- 
tion II and the experimental results presented in Sec- 
tion III that theory and experiment cannot be compared 
quantitatively. However, a number of qualitative com- 
parisons can be made pertaining to both high and low 
energy behavior. The transfer of spectral weight from 
high to low energies is clearly predicted for the hole doped 
systems [n < 1, Fig. 5(b) and (d) compared with Fig. 7(c) 
and (d)]. We note that this does not pertain to the two- 
magnon peak, which is not included in the theory, but 
rather to the systematic overall increase with hole doping 
over energies up to at least 1 eV for p > 0.05 (n < 0.95). 
On the electron-doped side for NCCO, the redistribu- 
tion is considerably weaker in the theoretical prediction 
[Fig. 5(a) and (c)]. In the experiment [Fig. 7(a) and (b)] 
there is little change at low energy similar to the theo- 
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FIG. 7: (Color online). High-energy Raman response of 
[(a),(b)] NCCO and [(c), (d)] LSCO for electron and hole dop- 
ing, respectively. The spectra are shown in [(a),(c)] Bi g and 
[(b), (d)] B2 g symmetry at various doping levels at tempera- 
tures of roughly 200 K. The peaks at low doping in (a) and (c) 
are due to two magnon scattering. In general with increasing 
doping level the spectra on the hole doped side lose intensity 
while the intensity of the spectra on the electron doped side 
increases. Part of the results are similar to those of other 
authors 51-53,55 or were already published in Ref. 49. 



retical prediction, but in opposition to theory an increase 
at higher energies which is most likely originating from 
luminescence as outlined in Sec. IIIB. A significant en- 
hancement of spectral weight is observed at low energies 
in Big symmetry as predicted in the theoretical results. 
The derivation of the exact energy dependence of 
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FIG. 8: (Color online). Temperature dependence of the low- 
energy Raman response of overdoped (a) NCCO and (b) 
LSCO in B2 g symmetry. Only for the electron doped com- 
pound does there develop a quasiparticle peak at low temper- 
atures. Results for NCCO with x = 0.15 displaying similar 
properties were published in Ref. 63. The data for LSCO are 
adapted from Ref. 49. 



the electron-hole continuum is experimentally challeng- 
ing since various processes may contribute. Among the 
intrinsic contributions are luminescence, resonance en- 
hancement, and spin excitations. There are various stud- 
ies at relatively low Raman shifts on the dependence of 
excitations on laser photon energy 6 whereas there is 
less material on the high-energy continua. Studies of 
Kang et al. 67 and Blumberg et al. 68 show that the in- 
fluence of resonances on the low-energy electronic part 
are mild for blue-green excitation. These are the wave- 
lengths used in the experimental portion of this study. 
In addition, resonances do not affect the form factors di- 
rectly since they are dictated by symmetry, but rather 
affect the relative intensities of the channels and their 
possible dependence on the excitation energy. However, 
a calculation of the full vertices is certainly beyond the 
scope of the present paper and we confine our argumen- 
tation to the symmetry part of the vertices. 

At low energy, as one can see in Fig. 6, the Hubbard 
model predicts relatively flat spectra on the hole-doped 
side [panels (b) and (d)] and well-defined peaks for elec- 
tron doping [panels (a) and (c)]. Since the temperature in 
the simulations is high, the peaks are wide; however, for 
the electron doped systems these peaks do originate from 
quasiparticle-like, intraband particle-hole excitations in 
an almost normal metallic band (see Fig. 3). Somewhat 
differently in the experiment, the peak at low energy ap- 
pears only at relatively low temperatures. At the mo- 
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ment we do not know what kind of interactions lead to 
a reduction of the carrier lifetime already around room 
temperature. Empirically, the possibility exists that the 
relatively large coefficient A in front of the T 2 term of 
the resisitivity (found at least for LSCO 69 ) leads to a 
rapid suppression of the coherence peak with increasing 
temperature. 

Why do we believe that this low-energy peak and the 
spectral weight suppression originate from particle-hole 
excitations and not from fluctuations and a pseudogap 
such as in Y-123 56 and LSCO 70 ' 71 ? Simulations 62 show 
that the low-energy peak and the dip can in fact result 
from a reduced quasiparticle damping in a Fermi liquid 
phenomenology. In addition, the initial slope of the spec- 
tra [Eq. (3)] follows the resistivity, at least qualitatively, 
as shown in Fig. 9(a) while the case is opposite for LSCO 
[see the B\ g response in Fig. 9(b)] and Y-123. 

As previously noted, the spectra of materials close 
to half filling corresponding to n ~ 1 or x ~ have 
prominent peaks from magnetic excitations which are 
not reproduced by the current theory as only the low- 
est order approximation, i.e. the bare bubble, has been 
used for the calculation of the Raman response. 32,33 Sys- 
tematic studies in Y-123, Bi 2 Sr 2 CaCu20 8+I 5 (Bi-2212), 
Tl 2 Ba 2 Cu0 6+( 5 (Tl-2201), and LSCO show 49 that this 
approximation is too simple, in particular at lower dop- 
ing. As demonstrated in Bi-2212 the collapse of the ap- 
proximation appears to set in rather abruptly at n ~ 0.79 
(p ~ 0.21) and has no direct correspondence in the single- 
particle spectral function. 60 

How might one expect the overall features to differ un- 
der a scenario in which the Hubbard U were allowed to 
decrease significantly with doping away from half-filling 
on either the hole- or electron-doped sides of the phase di- 
agram? One expects a more pronounced transfer of spec- 
tral weight from high to low energies than that observed 
in Fig. 5 together with a pronounced shift of residual high 
energy spectral weight to lower energies corresponding to 
the reduction in Hubbard U. Doping-dcpcndcnt U also 
would lead to the appearance of sharper, more metallic 
quasiparticle peaks at low energies especially on the hole- 
doped side of the phase diagram at greater variance with 
experimental observations compared to the current theo- 
retical analysis. While one can argue that vertex correc- 
tions may tend to reduce or flatten the Raman response 
at low energies and introduce features corresponding to 
magnetic excitations lacking in the lowest order approx- 
imation used in this study, the effect of these corrections 
with a doping-dependent U would tend to decrease sig- 
nificantly with doping and potentially suppress magnetic 
excitations in the response more rapidly than indicated 
by experimental observations. 



B. Relaxation Rates 

The method for extracting Raman relaxation rates 
from the spectra has been described elsewhere. 72 The 




Temperature (K) Temperature (t) 



FIG. 9: (Color online). Temperature dependence of the 
(a) and (b) experimental and (c) and (d) theoretical Ra- 
man resistivities in B\ g (red) and B2 g (purple) symme- 
try obtained from the electron (Nd2- x Ce x Cu04) and hole- 
doped (La2- x Sr x Cu04) samples around optimal doping, re- 
spectively. 



calculated Raman relaxation rates from the experimen- 
tal spectra are shown in Figs. 9 (a) and (b) in addition to 
the inverse initial slopes of the theoretical spectra as de- 
termined via Eq. (3) as shown in Figs. 9 (c) and (d). For 
both samples the experimental Raman relaxation rates 
are higher in the B\ g than in the _B 2ff channel. The 
relaxation rates in i? 29 have a similar temperature de- 
pendence for the electron as well as for the hole doped 
side showing a concave curve following approximately T a 
with 1 < a < 2. 

On the hole-doped side the relaxation rates in B\ g have 
a convex form (a < 1). This is different from the electron 
doped side. Here the form changes from convex to con- 
cave for a small variation in doping from x — 0.16 to 0.17. 
This is actually the range of doping where a crossover be- 
tween a small and a large Fermi surface is observed with 
quantum oscillations. 73 The inverse initial slopes of the 
theoretical spectra are different in that the relaxation 
rates in Bi g are bigger than in B\ g for all temperatures 
on both sides of the phase diagram. Additionally the 
temperature dependence is concave for the electron as 
well as for the hole doped side. However, what can be 
qualitatively resolved, in agreement with experiments, is 
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the fact that the relaxation rates for both symmetries are 
smaller on the electron than on the hole doped side of the 
phase diagram. 

The huge difference in the experimental relaxation 
rates between the B\ g and B^g response in LSCO is not 
yet fully understood. As mentioned above it appears in a 
doping range in which the dichotomy between nodal and 
anti-nodal quasiparticles is still small 15 while the Raman 
spectra exhibit an abrupt onset of strong relaxation close 
to (ir, 0). 60 The discrepancy appears to indicate that Ra- 
man vertices should be renormalized at least at low dop- 
ing. The observed theoretical behavior for hole-doping 
is in better qualitative agreement with the experimental 
results on overdoped La2- x Sr x Cu04 (not shown) where 
experimentally one observes more Fermi liquid-like be- 
havior with metallic quasiparticles in the normal state 
as opposed to the strange metal and pseudogap phases 
at doping levels closer to half-filling. At the accessible 
temperatures in the theoretical analysis, there is no indi- 
cation of a pseudogap in the single-particle spectra; how- 
ever, from other small cluster studies one anticipates the 
formation of a pesudogap at considerably lower tempera- 
tures. 74 The pseudogap manifests as a reduction in single- 
particle spectral weight near the (tt, 0) points and one 
anticipates a concomitant reduction in the B\ g Raman 
response and significant increase in the corresponding 
relaxation rate potentially bringing the theoretical and 
experimental results into better agreement. Lower tem- 
peratures also tend to sharpen quasiparticle features po- 
tentially reducing relaxation rates, especially in the B^g 
channel extracted from quasiparticles near ~ (7r/2,7r/2) 
that are particularly wide in the current simulations. 



C. Evolution of the spectral weight with doping 

Finally, we analyze the ratio of the Raman intensi- 
ties in the B\ g and i?2g channel /Big/^B2g- For deriving 
-feig/^B2g we integrate the experimental intensity in B\ g 
and B2g channels over the range 800-1000 cm -1 , labeled 
/eig and iB2g, respectively. Over this range the experi- 
mental spectra are fairly temperature independent. 72 For 
a simple tight-binding band structure, one expects this 
ratio to be given by (£/2i') 2 ; however, discrepancies be- 
tween the experimental response and this simple expec- 
tation, at least on the hole-doped side of the phase di- 
agram, were pointed out in earlier studies. 75,76 While a 
similar analysis of integrated intensity from the theoreti- 
cal results is complicated by the lack of vertex corrections 
and the relatively high simulation temperature leading to 
significant low energy peaks in the response rather than 
flat, featureless spectra in the indicated energy range (see 
Fig. 6), we perform a corresponding analysis of the inte- 
grated spectral weight as a function of doping and tem- 
perature over the same energy range (see Fig. 6) which 
at least removes a degree of arbitrariness from the com- 
parison between experiment and theory. 

Fig. 10 shows the ratio -?Big/-?B2g both calculated from 
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FIG. 10: Ratio of the intensities in /Bi g /^B2 g as extracted 
from the theoretical (upper) and experimental (lower) Raman 
spectra for electron and hole doping levels indicated in each 
panel. 



the theoretical response and extracted from the experi- 
mental spectra. On the electron-doped side the experi- 
mental ratio is slightly temperature and doping depen- 
dent for those samples analyzed near optimal doping. 
While the ratio increases slightly with decreasing doping, 
the temperature dependence of the ratio remains rela- 
tively unchanged, increasing moderately with increasing 
temperature. From the theoretical results, we analyze 
the ratio both for optimal doping as well as for under- 
doped and overdoped systems. Close to half-filling, we 
find that the ratio increases significantly with decreas- 
ing temperature although it remains within a factor ~ 2 
of the value extracted near optimal doping. While this 
temperature dependence differs significantly from that 
found experimentally near optimal doping, it is not com- 
pletely unexpected given the experimental and theoreti- 
cal evolution of the Fermi surface and bandstructure with 
doping. 10 At low doping levels small, electron-like Fermi 
pockets first appear near (±7r, 0) and (0, ±7r) at low tem- 
peratures, due to antiferromagnetic folding, which in the 
simple response calculated here would lead to a signif- 
icantly larger low energy response in the B\ g channel 
compared to the B^g channel. With increasing tempera- 
ture, one expects increasing response in the B-2 g channel 
and a reduced ratio due to an increase in spectral weight 
in the nodal region. With increased electron doping the 
spectral weight in the nodal region increases leading to a 
reduction in the intensity ratio near optimal doping. The 
theoretical ratio increases slightly with decreasing tem- 
perature, at odds with the experimental result; although 
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it appears to saturate at the lowest simulated temper- 
atures and has a value that compares quite well with 
that extracted from experiments. Upon additional elec- 
tron doping a well-defined hole-like Fermi surface forms, 
in agreement with experiment; however, the antinodal 
BZ crossing is pulled progressively away from the (ir, 0) 
points leading to a reduction in the integrated Raman re- 
sponse in the B\ g channel with a subsequent reduction in 
the intensity ratio. While the temperature dependence of 
the ratio extracted from the theoretical Raman response 
is not in quantitative agreement with the experimentally 
determined intensity ratio, the general reduction in the 
ratio with increased doping is captured by the theoretical 
analysis and one also should keep in mind the limitations 
of the theoretical analysis (lack of vertex corrections) and 
significant difference in temperature scales. 



On the hole-doped side we can compare data over a 
much wider doping range. For low doping, at the onset 
of superconductivity (x = 0.05), the experimental ratio 
increases slightly with increasing temperature, but has a 
value that remains less than one over the studied tem- 
perature range. Upon hole doping away from half-filling, 
one observes Fermi arcs experimentally in the single- 
particle spectral function centered on the nodal points 
with the appearance of a pseudogap in the antinodal re- 
gion. 15 The effect of these features would tend to sup- 
press the integrated B\ g Raman response and hence the 
intensity ratio within the simple picture used to describe 
the response theoretically, in apparent agreement with 
the experimental findings. With additional hole dop- 
ing, the Fermi arcs connect as the pseudogap closes in 
the single-particle spectral function near optimal doping. 
Experimentally, the integrated intensity ratio increases 
by a factor of ~ 2 — 3 over its value at underdoping and 
the ratio increases with decreasing temperature. Upon 
overdoping, the antinodal BZ crossings move toward the 
(tt, 0) points as the van Hove singularity approaches the 
Fermi level in the same region of the BZ. This would 
significantly increase the integrated B\ g response in the 
simple band picture used in the theoretical analysis. Ex- 
perimentally the intensity ratio for the overdoped system 
increases substantially with decreasing temperature and 
has a value ~ 3 times larger than that found near op- 
timal doping. The ratio extracted from the theoretical 
response agrees qualitatively with the experimental re- 
sults both in the general trends with hole doping away 
from half-filling and changes in the ratio with increas- 
ing or decreasing temperature. The lack of quantitative 
agreement may be attributed to the relatively high sim- 
ulation temperatures that preclude the appearance of a 
pseudogap in the single-particle spectral function at low 
doping and that significantly broaden quasiparticle fea- 
tures irrespective of the doping level and to the lack of 
vertex corrections that would renormalizc the Raman re- 
sponse in each channel, particularly at low doping. 



V. CONCLUSIONS 

While it is quite clear that the theoretical and experi- 
mental results presented here cannot be compared quan- 
titatively, qualitative comparisons exist that pertain to 
both the high and low energy Raman response in hole 
and electron doped systems. Speaking generally, the the- 
oretical response shows a significant transfer of spectral 
weight from high to low energies with increased hole dop- 
ing away from half- filling in agreement with experiment. 
With electron doping, this transfer is less pronounced, 
except in the underdoped regime, also in agreement with 
experiment. At lower energies, the model calculations 
predict relatively flat spectra on the hole-doped side of 
the phase diagram and well-defined peaks associated with 
quasiparticle-like excitations for electron-doped systems. 
The lack of vertex corrections in the theoretical response 
means that this general agreement does not apply to 
the two-magnon peak associated with magnetic excita- 
tions that appears prominently in the experiment, but 
rather to the overall, systematic trends in the data that 
appear with either electron or hole doping and already 
reveal an electron/hole doping asymmetry. While sim- 
ilar overall trends would be expected from a theoreti- 
cal analysis including significant reductions in electron 
correlations with doping away from half-filling (doping- 
dependent Hubbard U), the effect would be more pro- 
nounced and at greater variance with the experimental 
observations. In particular, one expects more metallic 
behavior at even lower electron or hole doping and would 
anticipate significantly reduced influence from vertex cor- 
rections with doping that would be necessary to bring 
the theoretical and experimental observations into better 
qualitative and quantitative agreement with one another. 

Additional electron/hole doping asymmetry has been 
observed experimentally in both the Raman relaxation 
rates in the B\ g and Bi g channels as well as in the ratio 
between the integrated Raman intensities in these two 
channels. The Raman relaxation rates or resistivities ex- 
tracted from the theoretical model near optimal doping 
reveal a similar asymmetry at least in the B\ g channel. 
As argued, the inclusion of pseudogap behavior in the 
antinodal region and sharper quasiparticlc-like features 
in the near nodal region with reduced temperatures may 
bring the theoretical and experimental results into better 
qualitative agreement, especially concerning the behav- 
ior of the Raman resistivity in the B^ g channel where 
there is a less pronounced asymmetry between electron 
and hole doping. 

The experimentally observed ratio between the inte- 
grated intensities in the B\ g and B^ g channel has been 
shown systematically over a wide range of hole doping 
and temperature and near optimal doping over a wide 
range of temperatures on the electron-doped side of the 
phase diagram. On the hole-doped side of the phase di- 
agram, the ratio extracted from the theoretical model 
agrees qualitatively with the experimental results both 
as a function of doping and temperature, although the 
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effective temperatures differ significantly. The relatively 
high simulation temperatures preclude the appearance 
of a pseudogap in the single-particle spectral function 
at low doping and broaden quasiparticle features, espe- 
cially at high doping. Vertex corrections, missing from 
the theoretical analysis, would significantly renormalizc 
the Raman response in each channel, particularly at low 
doping. Incorporating both of these corrections would 
improve the quantitative agreement between the theo- 
retical and experimental ratios in hole-doped systems. 

On the electron-doped side of the phase diagram near 
optimal doping the experimental intensity ratio is only 
slightly temperature and doping dependent. From the 
theoretical model we have analyzed the integrated in- 
tensity ratio over a wider range of electron doping and 
predict a significant difference between the electron and 
hole doping evolution. Close to half-filling, the theoret- 
ical ratio increases significantly with decreasing temper- 
ature, showing similar behavior to the theoretical and 
experimental ratio observed for hole-overdoping. With 
additional electron doping the ratio progressively de- 
creases. The doping evolution of the theoretical ratio 
can be partially understood by considering the evolution 
of the single-particle spectral function and the Fermi sur- 
face. 10 ' 15 While the Fermi surface consists of small Fermi 
arcs with a significant pseudogap on the hole-underpoded 
side of the phase diagram leading to a small integrated 
intensity ratio, the Fermi surface consists of small elec- 
tron pockets on the electron- undcrdoped side of the phase 
diagram leading to a rather large integrated intensity ra- 
tio. The evolution of the Fermi surface and bandstruc- 
ture with additional doping on cither side of the phase 
diagram appears to correlate quite well with the doping 
evolution and the electron/hole asymmetry observed in 
the intensity ratio. 

While the general trends observed in the theoretical 



and experimental results agree qualitatively and reveal 
significant electron/hole doping asymmetries in several 
quantities, one should keep in mind the limitations of 
the theoretical approach: (1) the simulation tempera- 
tures are quite high when compared to the experiments 
moaning that effects from features like the pseudogap 
that manifest at much lower energies are neglected, (2) 
the theoretical analysis based on the lowest order ap- 
proximation to the Raman response neglects important 
vertex corrections that renormalize the response, espe- 
cially at low energy, and therefore precludes the appear- 
ance of magnetic excitations that are prominent in the 
experiment. Although in general the results are encour- 
aging, further improvements including materials speci- 
ficity missing from the current theoretical analysis and 
the inclusion of vertex corrections by a direct determina- 
tion of the Raman response, typical for the charge and 
spin response functions (charge and spin dynamical struc- 
ture factors), would certainly improve agreement. Hence, 
apart from some qualifications explicitly stated, the Hub- 
bard model is found also here to capture the essential 
physics of the cuprates around half-filling. 
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